Structure and energetics of Si(lll)-(5x2)-Au 
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We propose a new structural model for the Si(lll)-(5x2)-Au reconstruction. The model incorpo- 
rates a new experimental value of 0.6 monolayer for the coverage of gold atoms, equivalent to six gold 
atoms per 5x2 cell. Five main theoretical results, obtained from first-principles total-energy calcu- 
lations, support the model. (1) In the presence of silicon adatoms the periodicity of the gold rows 
spontaneously doubles, in agreement with experiment. (2) The dependence of the surface energy on 
the adatom coverage indicates that a uniformly covered phase is unstable and will phase-separate 
into empty and covered regions, as observed experimentally. (3) Theoretical scanning tunneling 
microscopy images are in excellent agreement with experiment. (4) The calculated band structure 
is consistent with angle- resolved photoemission spectra; analysis of their correspondence allows the 
straightforward assignment of observed surface states to specific atoms. (5) The calculated activa- 
tion barrier for diffusion of silicon adatoms along the row direction is in excellent agreement with 
the experimentally measured barrier. 



I. INTRODUCTION 

Forty years ago Bishop and Riviere first observed the 
(111) surface of silicon to reconstruct, with fivefold pe- 
riodicity, in the presence of gold,^ Since that time the 
Si(lll)-(5x2)-Au reconstruction has been widely stud- 
ied, in several hundred publications, as a prototype lin- 
ear metallic chain system in which the physics of one- 
dimensional metals is approximately realizedi^ 

These investigations have provided very substan- 
tial insights into many aspects of Si(lll)-(5x2)- 
Aui^i^ Less successful have been the many attempts 
to use the clues provided by experiment to con- 
struct a complete structural model for this com- 
plicated reconstruction. Beginning with the early 

work of LeLay over a dozen different models have 
been proposed, 6. 7. 8. 9 . 10.11.12. 13.14. is. 16.17. is. 19.20. 21. 22. 23. 24 

All were eventually found to be inconsistent with the 
results of scanning tunneling microscopy (STM), angle- 
resolved photoemission spectroscopy (ARPES), or both. 

In this paper we propose a new structural model for 
Si(lll)-(5x2)-Au that is fully consistent with all exper- 
imental data to which we have compared. The model 
is similar to one proposed by Erwin in 2003 but is 
modified to be consistent with a recently revised value of 
the gold coverage.^- The modifications, although seem- 
ingly minor, for the first time bring the predictions 
of the model— for STM and ARPES as well as other 
phenomena — into excellent agreement with experiment. 
More importantly, the new model opens the door to 
a more fundamental physical understanding of Si(lll)- 
(5x2)- All based on its detailed atomic structure. 



II. STRUCTURAL MODEL 

The model proposed here is shown in Fig. [TJ The basic 
structure is similar to one proposed several years ago in 



Ref. '2^ but there are several important differences. For 
this reason it is useful to discuss both the similarities and 
differences between the older model and the one proposed 
here (hereafter the "2003 model" and the "2009 model"). 

The starting point for the 2003 model was the experi- 
mental observation that 0.4 monolayer (ML) of gold in- 
duces a stable reconstruction of the Si(lll) surface. Us- 
ing this coverage value the 2003 model was constructed 
with fo ur g old atoms per 5x2 cell. As shown in Fig. 1 
of Ref. [20, the gold atoms substitute for silicon atoms 
in the topmost surface layer, forming two Au-Si chains 
oriented along the [110] direction. Adjacent to these two 
chains is a silicon "honeycomb chain," a thin graphitic 
strip of silicon that owes its stability to a Si=Si double 
bond. The basic reconstruction just described has 5x1 
periodicity. The experimentally observed 5x2 substrate 
periodicity was argued to arise from a row of silicon re- 
bonding atoms that bridges a channel between one of the 
Au-Si rows and the silicon honeycomb chain. A full row 
of rebonding atoms overcoordinates some atoms, while 
a half-occupied row, with 5x2 periodicity, leaves some 
dangling bonds unsaturated. The latter arrangement was 
argued to be energetically preferred when the system is 
doped with extra electrons. These were supplied by sili- 
con adatoms adsorbed on top of the Au-Si chains. 

The recently revised experimental determination of the 
gold coverage as 0.6 ML obviously calls for a revised 
structural model as welli^ The 2009 model accommo- 
dates the additional two gold atoms per 5x2 cell by re- 
placing the half-occupied rebonding row of silicon atoms 
from the 2003 model with a full row of gold atoms. The 
other two building blocks of the 2003 model — the silicon 
honeycomb chain and the adsorbed silicon adatoms — are 
unchanged in the new model. The total coverage of sil- 
icon atoms in the top layer of the 2009 model varies, 
according to the coverage of adatoms, between 1.20 ML 
(for the undecorated surface) and 1.25 ML (for saturation 
coverage). These values are consistent with the range 
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of experimentally determined silicon coverage, 1.1 to 1.3 
ML (Ref.i26i), 1.23±0.003 ML (Ref.i221), and 1.3±0.1 ML 
(Ref. [H). The 2009 model is energetically more favor- 
able than the 2003 model: the theoretical surface energy 
is 17.1 meV/A^ lower in the Au-rich limit, in which the 
chemical potential is taken as the bulk energy per atom.^^ 

Much of the rest of this paper will focus on the inter- 
esting role played by the new row of gold atoms. But 
it is worth pausing briefly to place Si(lll)-(5x2)-Au in 
the larger context of other metal-induced reconstructions 
of Si(lll). It is now known that a great variety of 
metal adsorbates induce closely related reconstructions 
based on the silicon honeycomb chain. But the details 
differ, sometimes with unexpected consequences. Eas- 
ily ionized adsorbates such as alkalis, alkaline earths, 
and some rare earths form a family of "honeycomb- 
chain channel" (HCC) reconstructions in which the ad- 
sorbates occupy channels between adjacent silicon honey- 
comb chainsji22i2iiH Each adsorbate is three-fold coordi- 
nated by silicon atoms in the honeycomb chain, which are 
too far away (3.0 A) to form covalent bonds, consistent 
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FIG. 1: (color online). Proposed structure of Si(lll)-(5x2)- 
Au with gold coverage equal to 0.6 monolayer. Large yellow 
circles are gold, small circles are silicon. The surface layer 
consists of a gold single row (S), gold double row (D), and sil- 
icon honeycomb chain (HC). The surface energy is minimized 
when this surface is decorated by silicon adatoms (dark blue) 
with 5x4 periodicity, as shown. In the presence of adatoms 
the 5x1 periodicity of the underlying substrate spontaneously 
doubles to 5x2 (gray outline) due to dimerization within the 
gold double row. 



with a picture of ionic charge donation and electrostatic 
attraction. The interactions between adsorbates are also 
electrostatic, but repulsive.— Within this family of HCC 
reconstructions the adsorbate coverage is determined by 
a simple electron-counting rule first proposed by Lee et 
alM-t^ and later generalized by Battaglia et alM--^ 

Adsorbates that are less ionic, such as silver and gold, 
also form reconstructions based on the silicon honey- 
comb chain but the role of the adsorbate is more interest- 
ing. Silver adsorbates occupy the channel of a HCC-like 
reconstruction,'^'^ but the stronger interaction between 
silver and silicon leads to a preference for two-fold coordi- 
nation of silver by silicon, at a much reduced distance of 
2.6 A. The two-fold coordination brings silver adsorbates 
sufhciently close to each other to allow pairing into sil- 
ver dimers. Chuang et al.^^ and Urbieta et al.^^ showed 
that the phase of this pairing alternates between adja- 
cent channels, modulating the periodicity of the basic 
3x1 HCC reconstruction to a c(12 x 2) variant. 

Returning now to Si(lll)-(5x2)-Au, it appears that 
the role of the adsorbate is still more complex. Gold is 
very reactive on silicon surfaces 4£ This reactivity is al- 
ready evident in the Au-Si rows of Fig. [1] Within these 
rows gold completely substitutes for the top layer of the 
surface silicon bilayer, with each gold atom covalently 
bonded to three silicon atoms. Likewise, each silicon 
atom at the edge of the silicon honeycomb chain forms a 
covalent bond to a gold atom. The Au-Si bond lengths, 
both within the Au-Si rows and bonded to the silicon 
honeycomb chain, are 2.5 A, smaller than for any other 
adsorbate studied. 

The most interesting aspect of the 2009 model is the 
behavior of the gold double row labeled "Au D" in Fig. [TJ 
As shown in the figure, the equilibrium geometry of this 
double row is dimerized. Period doubling was also part 
of the 2003 model, but its origin — the half-occupied re- 
bonding row — was simpler. In the 2009 model the dimer- 
ization occurs only in the presence of silicon adatoms, or 
when the surface is doped with extra electrons. In Sec. lIVI 
we show that these two scenarios are largely equivalent. 
We further demonstrate that in the presence of adatoms 
the dimerization is driven by an unusual "double" Peierls 
mechanism, in which the distortion opens two gaps in the 
band structure of the undistorted 5x1 substrate. 



III. METHODS 
A. Theoretical methods 

First-principles total-energy calculations were used to 
determine equilibrium geometries and relative energies 
of the basic model and its variants. The calculations 
were performed in a slab geometry with four layers of 
Si plus the reconstructed top surface layer and a vac- 
uum region of 8 A. All atomic positions were relaxed, 
except the bottom Si layer and its passivating hydrogen 
layer, until the largest force component on every atom 
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was below 0.01 eV/A. Total energies and forces were cal- 
culated within the PBE generalized-gradient approxima- 
tion to density-functional theory (DFT) using projector- 
augmented-wave potentials, as implemented in vASP4iii^ 
The plane-wave cutoff for all calculations was 250 eV. 

The sampling of the surface Brillouin zone was chosen 
according to the size of the surface unit cell and the rel- 
evant precision requirements. For example, the depen- 
dence of the total energy on dimerization (Fig. [2]) was 
calculated using a 5x4 unit cell and 2x2 zone sampling, 
with convergence checks using 4x4 sampling. The de- 
pendence of the relative surface energy on silicon adatom 
coverage (Fig. |3]) requires greater precision because the 
energy variations are smaller. Hence these surface ener- 
gies were calculated using a 5x8 unit cell (to allow an 
adatom coverage of 1/8) and 8x4 sampling, with conver- 
gence checks using 12x6 sampling. Finally, the potential- 
energy surface for adatom diffusion (Fig. [7|) was calcu- 
lated using a 5x4 unit cell and 2x2 zone sampling. 

Simulated STM images (Fig. |4]) were calculated using 
the method of Tersoff and Hamann.^'^ For the filled-state 
image we integrated the local density of states (LDOS) 
over a chosen energy window of occupied states up to the 
Fermi level; for the empty-state image the integration 
was over unoccupied states starting at the Fermi level. 
The simulated STM topography under constant-current 
conditions was obtained by plotting the height at which 
the integrated LDOS is constant. 
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FIG. 2: (color online). Variation of the total energy as a 
function of dimerization along the chain direction. The di- 
mensionless dimerization parameter is d = (ai— ao)/ao, where 
a\ is shown in the inset and ao is the surface lattice constant. 
Total energies were calculated with full relaxation for each 
constrained value of d. 



extra electrons are present. Next we examine how the 
surface energy varies as a function of silicon adatom cov- 
erage; we find that the energy is minimized for the cover- 
age 1 /4 shown in Fig. [T] Finally, we show that the model 
explains the existence of a localized defect structure seen 
even on carefully prepared surfaces. 



B. Experimental methods 

Silicon wafers (from Virginia Semiconductors) were de- 
gassed for several hours at 700 °C before flashing at 1250 
°C for a few seconds. A rapid cool-down to 850 °C was 
followed by slow cooling to room temperature. An im- 
portant prerequisite for obtaining well-defined row struc- 
tures was a flexible mount of the samples that prevented 
strain from building up during high-temperature flashing. 
Gold was evaporated from a Mo wire basket. The pres- 
sure was kept below 5x 10^" mbar throughout the sample 
preparation. All STM measurements were carried out at 
room temperature with low tunneling currents (<50 pA). 

Band dispersions were obtained using a Scienta 200 
spectrometer with E', d mult idetect ion and an energy res- 
olution of 20 meV for electrons and 7 meV for photons. 
We used p-polarized synchrotron radiation at a photon 
energy hv — 34 eV, where the cross section of silicon sur- 
face states has a maximum relative to the bulk states. 
The data are for a sample temperature below 100 K. 



IV. ENERGETICS 

This section addresses three issues related to the en- 
ergetics of Si(lll)-(5x2)-Au. We consider first the ener- 
getics of dimerization, and show that dimerization occurs 
naturally within the 2009 model when silicon adatoms or 



A. Dimerization of the substrate 

The model shown in Fig. [T] has 5x4 periodicity, be- 
cause silicon adatoms decorate the surface in a 5x4 ar- 
rangement. But it is clear from the figure that the un- 
derlying substrate (that is, ignoring the adatoms) can be 
understood more simply as a 5x1 reconstruction whose 
periodicity is doubled, along the chain direction, to 5x2 
by dimerization within the Au double row. Understand- 
ing the nature and origin of this dimerization is impor- 
tant for explaining many experimental aspects of Si(lll)- 
(5x2)-Au, including the fine details of STM imagery, 
data from ARPES, the observation of nanoscale phase 
separatiouf^^i^i^i^i^ the diffusion of silicon adatoms on 
the surfacCji ^^ ii " " ' and the existence and motion of domain 
walls within the Au-Si rowsi^ 

We begin by defining more precisely the nature of the 
dimerization. The inset to Fig. |2| shows a detail of the 
gold double row near a silicon adatom. Within this 
double row the gold atoms have a ladder-like arrange- 
ment, with Au-Au bonds as the rungs. In the absence of 
adatoms (or extra electrons) these rungs are all parallel 
to each other, with a spacing equal to the silicon surface 
lattice constant ao. When adatoms decorate the surface 
in the 5x4 arrangement shown in Fig. [1] the rungs rotate 
away from their parallel alignment. This rotation occurs 
almost completely within the (111) surface plane, and 



the rungs are quite rigid: the Au-Au bond length (2.94 
A) changes by less than 1%. This rigidity is not sur- 
prising, because the Au-Au bond length is already very 
close to the bulk gold bond length (2.88 A). The sign of 
the rotation alternates along the chain direction. Hence 
the dimerization can be viewed as an antiferrodistortive 
instability. 

To quantify the dimerization one could, of course, use 
the angle of rotation of the Au-Au rungs. We choose 
instead a more physically transparent measure: the dis- 
tance ai between gold atoms on the side of the double 
row adjacent to the silicon honeycomb chain, as labeled in 
Fig. O A dimensionless dimerization parameter can then 
be defined, d = (ai — ao)/ao. For the 5x4 arrangement 
of adatoms shown in Fig. [l] the equilibrium dimerization 
parameter is dcq — 0.14. 

Some insight into the origin of the dimerization may 
be obtained by computing the DFT total energy Et as 
a function of d while relaxing all other degrees of free- 
dom. The results are shown in Fig. [2] for three variants 
of the basic model: the undecorated and undoped surface 
(without adatoms or extra electrons); the adatom-doped 
5x4 surface of Fig.[Tl and the undecorated surface doped 
with two extra electrons per 5x4 cell. For the undeco- 
rated undoped surface the minimum is at d = 0, indi- 
cating that dimerization is not stable. For the adatom- 
doped surface there is a single minimum in the energy, as 
expected, at d — +0.14. For the electron-doped surface 
the behavior of Et{d) is nearly indistinguishable, for posi- 
tive d, from that of the adatom-doped surface. This sim- 
ilarity strongly suggests that each silicon adatom dopes 
two electrons to surface states. More substantive evi- 
dence for this conjecture is found in the electronic struc- 
ture of these variants, as we show below in Sec. IVII 

The behavior of Et{d) for negative d, however, is very 
different. This is because adatom-doping strongly breaks 
the 5x2 symmetry of the surface, while electron-doping 
does not. This suggests another role for the silicon 
adatoms: to pin the phase of the dimerization such that 
d > at the position of the adatom. We will return to 
this role in Sec. IVIII when considering the diffusion of 
adatoms. 



B. Adatom coverage and phase separation 

The preceding discussion and calculations were based 
on the assumption that silicon adatoms decorate the 
surface in a 5x4 arrangement, equivalent to a cover- 
age of 1/4 adatom per 5x1 cell. This phase can be 
achieved experimentally by depositing silicon onto the 
surface at temperatures around 300 °C until saturation 
is reached. But this saturated phase is metastable: 
annealing causes half of the silicon adatoms to diffuse 
away. The resulting equilibrium phase is not uniform, 
but instead exhibits patches with local adatom cover- 
age of 1 /4 interspersed with patches of undecorated sur- 
face; the global average adatom coverage is close to 
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Si adatom coverage (per 5x1 cell) 

FIG. 3: Theoretical surface energy as a function of silicon 
adatom coverage. Energies were calculated at the labeled cov- 
erages; the interpolating curve is a guide to the eye. Adatoms 
occupy their preferred binding sites as shown in Fig. [1] The 
surface energy is lowest for coverage 1/4, corresponding to 
the fully saturated 5x4 arrangement of adatoms. The dot- 
ted line highlights a local maximum of the energy, and shows 
that a surface with coverage less than 1/4 will phase separate 
into a mixture of empty and 1/4-covered regions, as found 
experimentally. 

-[^^g^44j45^j47j48 ^j^jg experimental result poses two ques- 
tions for theory. (1) What coverage of silicon adatoms 
minimizes the calculated surface energy? (2) Does the 
observed patchiness of the adatom distribution arise from 
an instability toward phase separation? 

To compare the energies of phases with different silicon 
adatom coverage, we compute the relative surface energy 

Es ^ EtiNsi) ~ Nsi fxsu (1) 

where Et{Nsi) is the total energy of a surface unit cell 
containing Nsi silicon atoms (including adatoms), and 
/isi is the silicon chemical potential. Thermodynamic 
equilibrium between the surface and the bulk requires 
/iSi to be the energy per atom in bulk silicon. Although 
we do not explicitly consider the free energy at finite 
temperature, configurational entropy effects may indeed 
play a role; see below. 

Within this formalism, adatoms are thermodynami- 
cally stable only if their presence lowers the surface en- 
ergy of the undecorated surface. Equation [T] shows that 
this requirement implies that the adsorption energy per 
adatom must be greater than the silicon cohesive energy 
(5.4 eV within DFT/PBE). Adsorption energies below 
this threshold imply that adatoms may be temporarily 
metastable, but not thermodynamically stable. This is a 
very stringent requirement, and one that is rarely satis- 
fied in our experience (it was not in the 2003 model). 

We calculated the relative surface energies for five dif- 
ferent silicon adatom coverages. To minimize numerical 
uncertainties the same 5x8 supercell was used for each 
coverage, and full relaxation performed in every case. 
The results are shown in Fig. [3] The surface energy 
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FIG. 4: Comparison of experimental and simulated (inset) STM images for Si(lll)-(5x4)-Au. (a) Filled states, V = —0.7 eV. 
(b) Empty states, V = +1.0 eV. The simulated images were obtained using slightly different voltages: V = — f.O and +0.5 
eV, respectively. Projected positions are shown for silicon adatoms and top-layer silicon and gold atoms; the size of the circles 
indicates the height of the atoms. 



is minimized for a coverage of 1/4. This is the cover- 
age depicted in Fig. [U and agrees well with the experi- 
mentally observed local coverage within adatom-covered 
patches. In Sec. lVII we show why the value 1/4 is special: 
at this adatom coverage — but at no other — the surface 
band structure is fully gapped and hence favored because 
occupied states can move down in energy. 

We turn now to the second question posed above. 
Of the five coverages considered here, the second most- 
favorable is not the adjacent 1/8 or 3/8 phase, but rather 
the undecorated "empty" surface. In other words, the 
surface energy of the intermediate 1/8 phase is higher 
than the average of the two endpoint phases, empty and 
1/4. This result suggests that a range of such intermedi- 
ate coverages between and 1 /4 may have energies above 
the tie line, shown in Fig. [31 connecting the two endpoint 
phases. This implies that a surface prepared with adatom 
coverage 1/8 (and perhaps any intermediate coverage be- 
tween and 1/4) will phase separate into a mixture of 
empty and l/4--covered regions. 

Of course, this conclusion does not answer the related 
question of why the observed global average coverage is 
close to 1/8. We hypothesize that the 1/8 phase may 
be stabilized at finite temperature by the entropy gained 
from occupying only half the adatom sites of the 1/4 
phase. We also do not address here the characteristic 
size of the phase-separated regions; this would require 
analyzing the energy cost of forming the phase bound- 
ary between the and 1/4 phases. Finally, we leave for 
future analysis the possibility that the phase-separated 
state benefits electrostatically from the charge transfer 
proposed by Yoon et al. to take place between the and 
1/4 phases4i 



C. Low-energy defect structure 

In addition to the phase separation just discussed, 
STM images of Si(lll)-(5x2)-Au sometimes reveal the 
presence of an occasional "mirror-domain" defect.— 
The signature of this defect is a slight shift, by roughly 
3 A to the right in Fig. (TJ of the bright protrusions that 
dominate the STM topography. The intensity of the pro- 
trusions and the position and shape of the underlying 
features are not affected. 

A simple variation of the basic structural model ac- 
counts for the observed defects. Starting from the struc- 
ture shown in Fig. (TJ one can construct a mirror-image 
variant by reflecting the entire top layer of atoms through 
the vertical plane bisecting the Si=Si double bonds of the 
honeycomb chain. (This mirror plane was discussed in 
Ref. [S^ and is shown in Fig. 1 therein.) The reflection 
leaves the honeycomb chain unaffected but reverses the 
ordering and orientation of the Au-Si rows. The geome- 
try of the rigidly reflected structure is already extremely 
close (within a few hundreths of an Angstrom) to the 
fully relaxed geometry. The reflection shifts the silicon 
adatoms by 2.2 A to the right, consistent with the ob- 
served shift. 

The surface energy of the relaxed defect reconstruc- 
tion is only slightly higher, by 1.2 meV/A^, than that 
of the original reconstruction. From inspection of the 
atomic positions of the original and defect models, we 
judge the energy cost of forming an interface between 
the two phases to be quite small. Thus it is plausible 
that short sequences of this defect structure — manifested 
as a small rightward shift of the silicon adatoms — should 
appear even on carefully prepared surfaces, with mini- 
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mal disruption to either the primary reconstruction or 
the local adatom coverage. 



V. COMPARISON WITH STM 

Two decades of studies probing the Si(lll)-(5x2)- 
Au surface with STM have created a richly detailed 
and consistent picture of its topography and bias 
dependence.— '^'^i^^'^^i^'^'^^'^^i^^ These data make possible 
a very stringent test for any proposed structural model 
by comparing theoretically simulated STM images to 
atomic-resolution experimental images. In this section 
we show that simulated images based on the 2009 model 
are in excellent and detailed agreement with recent ex- 
perimental images. Moreover, the new model resolves 
a small but irritating discrepancy found in comparisons 
based on the earlier 2003 model. 

The most pronounced features in STM imagery of the 
Si(lll)-(5x2)-Au surface are the bright protrusions now 
well established as originating from the silicon adatoms. 
In topographic maps these are generally the highest 
features in both filled- and empty-state images. As 
discussed elsewhere**^ and in the previous section, the 
adatoms occupy a 5x4 lattice. Much detailed informa- 
tion is contained in the imagery in between these lattice 
sites; the topography of this substrate has 5x2 periodic- 
ity. Of special interest is the registry of the 5x4 adatom 
lattice and the 5x2 substrate; the relative alignment of 
these two lattices provides an important test for struc- 
tural models. (We do not address here the lack of corre- 
lation between the adatoms in different rows, previously 
discussed in Ref. lisl .) 

Figure m shows experimental and theoretical simulated 
STM images for filled and empty states of the Si(lll)- 
(5x2)-Au surface, in a region where the silicon adatom 
coverage is 1/4. The agreement between experiment and 
theory is excellent, and allows all of the main experimen- 
tal features to be easily identified. (1) The dark ver- 
tical channels separating the three main rows shown in 
Fig. [4] arise from the doubly-bonded rungs of the silicon 
honeycomb chains. (2) The bright protrusions are from 
the silicon adatoms. (3) The triangular features with 
x2 periodicity, constituting the right edge of the rows, 
arise from the combination of two contributions: a pair 
of outer gold atoms brought close by dimerization (the 
apex of the triangle) and two silicon atoms at the left 
edge of the honeycomb chain (the base of the triangle). 
(4) In the filled states, the left edge of the rows is not well 
resolved, and appears to have either xl or very weak x2 
periodicity; these spots arise from the silicon atoms at 
the right edge of the honeycomb chain. (5) In the empty 
states, the x2 triangular features on the right edge of 
the rows alternate with V-shaped features that open to 
the left. These are a combination of two contributions: 
a pair of dimerized inner gold atoms (the apex of the 
V) and Si-Au bonds that are slightly dimerized by their 
proximity to the gold double row (the arms of the V). 



In the experimental images the alignment of the 5x4 
adatom lattice and the 5x2 substrate topography is as 
follows. The bright protrusions are located slightly to 
the left side within the main rows — except in the defect 
regions discussed in the previous section, where the pro- 
trusions are shifted to an equivalent location on the right 
side of the row. This off-center location is accurately 
reproduced in the simulated images based on both the 
original model and the defect model (not shown). Along 
the rows of the experimental images, the 5x4 protrusions 
are symmetrically straddled by the 5x2 triangular and 
V-shaped features of the substrate. Fig.|3]shows that this 
symmetric registry is properly reproduced in the simu- 
lated images — resolving a problem with the 2003 model 
first pointed out by Seifert whereby the alignment of the 
two lattices was asymmetrici^l 



VI. ELECTRONIC STRUCTURE 

Photoemission studies of Si(lll)-(5x2)- 
Au began in the mid-199Gs and continue 
today^^i^i^i^i^i^i^i^'^iSOiei xhese studies have 

led to important insights into the electronic structure 
of this complicated surface and provide tests comple- 
mentary to STM for evaluating structural models. In 
this section we address three aspects of the electronic 
structure of Si(lll)-(5x2)-Au: the detailed mechanism 
that drives the period doubling discussed in Sec. IIVI 
a comparison of the theoretical band structure to 
angle-resolved photoemission data; and an explanation 
of which specific surface orbitals give rise to the observed 
band structure. 



A. Origin of substrate period doubling 

In earlier sections the geometry and energetics of the 
period doubling was explored without any consideration 
of the underlying mechanism. Here we suggest an ex- 
planation, based on the electronic structure, for why the 
5x1 substrate dimerizes to 5 x 2 in the presence of silicon 
adatoms or extra electrons. 

We begin by considering a variant of the full 5x4 model 
of Fig. [T]in which the silicon adatoms are absent. When 
this surface is relaxed within DFT the dimerization van- 
ishes and hence the periodicity reverts to 5x1. The the- 
oretical band structure for this surface is shown in Fig. 
[SJa) for wave vectors parallel to the chain direction and 
energies near the projected band gap. The bands are 
plotted in the Brillouin zone of the full 5x4 model so 
that comparison with the bands of the full model can 
later be made. The folding of the 5x1 bands into the 
5x4 zone creates degeneracies at the zone center F and 
the 5x4 zone boundary A4; the most important zone- 
boundary degeneracy is marked by an arrow. The folding 
also creates band crossings, the most important of which 
(circled) is inside the band gap, about one-fourth of the 
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FIG. 5: Calculated band structure of Si(lll)-Au in four different scenarios, (a) 5x1 undecorated undoped surface, (b) 5x2 
electron-doped surface (2 electrons per 5x4 cell), (c) 5x4 adatom-doped surface (1 adatom per 5x4 cell); (d) 5x4 adatom- 
doped surface with spin-orbit coupling included. In each panel the bands are plotted in the same 5x4 surface Brillouin zone. 
The Fermi levels in panels (b), (c), and (d) are set to zero, and the zone-boundary degeneracies in (a) and (b) are aligned in 
order to highlight the evolution of the bands with doping. Projected bulk silicon bands are shown in gray. The formation of 
a hybridization gap at a band crossing (circled) and the opening of a gap at the zone boundary (arrows) are indicated. The 
observed surface is phase-separated into undecorated 5x2 regions and adatom-doped 5x4 regions in equal proportion. 



way from F to A4. The Fermi level is very close to the 
top of the valence band and the system is metallic. 

Next we consider how this band structure changes 
when we dope the surface with extra electrons. It was 
demonstrated in Fig. [5] that a doping level of 2 electrons 
per 5x4 cell leads to stable 5x2 dimerization, with dis- 
tortion energetics nearly identical to that from silicon 
adatoms at 1 /4 coverage. Figure Efb) shows the band 
structure from this electron-doped surface. The electron 
doping has two important consequences: the antiferrodis- 
tortive dimerization creates a large hybridization gap at 
the band crossing of the the undoped surface, and the 
Fermi level is pushed into this gap. The band degener- 
acy at the 5x4 zone boundary remains intact (see arrow), 
because the electron-doped surface has perfect 5x2 peri- 
odicity. 

In the presence of silicon adatoms at 1 /4 coverage this 
last degeneracy is lifted. The last two panels of Fig. O 
show the bands calculated at two levels of theory: (c) 
scalar relativistic; (d) fully relativistic with spin-orbit 
coupling. (We will show in the next section that the 
gold character of these surface bands is substantial, hence 
the spin-orbit splitting is large, about 0.2 eV.) Panel (d) 
shows that the system now develops a full gap. At 1/4 
adatom coverage the Fermi level falls just inside this gap, 
making the system insulating. 

To summarize, we find that silicon adatoms at cover- 
age 1 /4 create a multiband metal-insulator transition on 
Si(lll)-(5x2)-Au. The first (electronically induced) gap 
arises from band hybridization originating from dimeriza- 
tion and the resultant lowering of symmetry. The second 
(adatom induced) gap arises from the 5x4 potential of 
the adatoms, which lifts the degeneracy at the 5x4 zone 
boundary A4. 



B. Comparison with angle-resolved photoemission 

The above discussion focused on bands above the 
Fermi level. Now we turn to the occupied states, where 
we can make direct comparison to experimental data. 

The experimental band structure is highly one- 
dimensional near the Fermi level, and becomes gradually 
more two-dimensional for lower energies; here we limit 
our discussion to the one-dimensional dispersion along 
the chain direction. The results of our ARPES studies of 
Si(lll)-(5x2)-Au are summarized in Fig. [Sl^a), which is 
reproduced from an earlier publication."*^ By combining 
momentum- and energy-distribution curves from several 
Brillouin zones, three bands can be identified: 1', 1", 
and 2. (These labels are used by analogy to comparable 
bands at stepped Si(lll) surfacesi^) Figure [Hfa) shows 
the bands in the repeated-zone scheme of a surface with 
5x1 periodicity, for which Ai is the zone boundary. The 
5x2 period doubling introduces backfolded replicas of 
these bands, shown as dashed lines. The backfolded band 
1' is indeed found where expected, indicating that this 
band has strong 5x2 character. The backfolded replica of 
band 2 is too weak to be observed, indicated that it has 
mainly 5x1 character. Band 1" becomes more intense as 
the silicon adatom coverage is increased, indicating that 
it has strong 5x4 character 42. 

It is difficult to compare directly the theoretical bands 
of the adatom-doped surface to these experimental re- 
sults, because the 5x2 bands of Fig.lHI^a) must be folded 
once more into the Brillouin zone of the 5x4 cell. This 
folding creates many bands in a small energy interval and 
overly complicates the comparison between theory and 
experiment. We choose instead a simpler and clearer ap- 
proximate approach: to compare the experimental 5x2 
bands of Fig. [DJa) to the theoretical 5x2 bands of the 
electron-doped surface with no silicon adatoms. In doing 
so it must be kept in mind that the 5x4 potential of the 
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FIG. 6: (color online) Comparison of angle-resolved pho- 
toemission data to the theoretical band structure. (a) 
Photoemission-derived band dispersion reproduced from Ref. 
|4^ . Band l' (red) has strong 5x2 periodicity, band 2 (yel- 
low) has mainly 5x1 periodicity, and band l" (gray) has 5x4 
periodicity and becomes more intense with increasing silicon 
adatom coverage.— Inset: Brillouin zones for 5x1, 5x2, and 
5x4 surface unit cells, (b) Theoretical band structure of the 
5x2 electron-doped surface. The diameter of each circle is 
proportional to the contribution from surface atoms in the 
Au-Si chains of Fig. [1] (c) Theoretical band structure of the 
hypothetical undecorated undoped 5x1 surface; these are to 
be compared with the folded bands in panels (a) and (b) . The 
labeled surface bands originate from gold and silicon orbitals 
in the single (S) and double (D) rows marked in Fig. 1. The 
double row leads to two bands, consisting of bonding (Db) 
and antibonding (Da) combinations of orbitals. 



adatoms, already seen to play an important role for the 
unoccupied bands, will be absent. 

Figure [6fb) shows the calculated bands for the 5x2 
electron-doped surface, plotted in the 5x1 repeated-zone 
scheme used in panel (a). (The projected bulk bands, 
however, are shown for reasons of clarity in the extended- 



zone scheme.) The diameter of each circle is proportional 
to the summed projections of the state onto gold and sil- 
icon atoms in the single and double Au-Si chains of Fig. 
[T] The solid colored curves represent our best effort to 
match the three strongest surface bands to the ARPES 
bands. The overall agreement for bands 1' and 2 is excel- 
lent, despite the complexity of the calculated bands even 
for this simplified surface. Note that in our interpreta- 
tion, 1' and 2 are not simple parabolic bands as depicted 
in Fig.[6l[a). Instead, each comprises two or more bands 
and exhibits several avoided crossings. Further support 
for this interpretation will be presented in the next sub- 
section, where we discuss the orbital origin of the bands. 

The agreement appears less satisfactory for band 1". 
Although this band is correctly centered at the A2 point, 
and the shallow dispersion near that point reasonable, it 
is shifted rigidly up in energy by 0.2 eV compared to ex- 
periment. We believe that this shift is a spurious effect 
arising from the omission of silicon adatoms in the cal- 
culation. Indeed, Choi et al. have shown experimentally 
that the 1" band shifts upward as the adatom coverage 
is reduced."*^ For the range of coverages studied (between 
37 and 97% of the saturation 1/4 coverage) the shift in 
energy was linear. Extrapolating this result to a surface 
free of adatoms, one naively expects an upward shift in 
energy of the 1" band by 0.21 eV. Such a shift would 
bring the results of ARPES and theory into excellent 
agreement for all three bands. 



C. Orbital origin of the bands 

We now make one last theoretical simplification by 
eliminating the extra electrons. Upon relaxation the un- 
doped surface is no longer dimerized, and the period- 
icity reverts to 5x1. This simplification is justified if 
the perturbation of the bands from the dimerization and 
the extra electrons is small. In the limiting case of zero 
dimerization the exact 5x2 bands can be obtained from 
the 5x1 bands by simple zone folding. Figure [SJc) shows 
the calculated band structure for the simplified surface, 
plotted in the conventional 5x1 reduced-zone scheme. It 
is readily apparent that the bands in (b) can indeed be 
accurately obtained from those in (c) by first folding the 
5x1 bands about the zone midpoint A2 , and then shift- 
ing the Fermi level upward by the appropriate amount 
(0.3 eV) for electron doping of rigid bands. 

The advantage of the simplified 5x1 band structure 
is that its orbital origin is easy to analyze, because the 
bands fall almost entirely in the projected gap and rarely 
cross. By examining individual states in real space, we 
find that the middle band S originates primarily from 
gold and silicon orbitals in the single chain "Au S" in Fig. 
[H The other two bands are more complicated. They arise 
from the double chain labeled "Au D." The two bands are 
the low-lying bonding (Db) and higher- lying antibonding 
(Da) combinations of orbitals of the gold atoms that con- 
stitute each rung of the ladder (combined with orbitals 
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FIG. 7; (color online). Theoretical potential energy surfaces 
for the diffusion of a silicon adatom along the chain direc- 
tion. Two approximations were considered for determining 
the diffusion pathway. (1) Substrate atoms were fixed rigidly 
at their equilibrium positions attained when adatoms are in 
their lowest energy site (circles). (2) Substrate atoms were 
allowed to conform to their "instantaneous" equilibrium po- 
sitions as the adatom diffused (squares). These two limiting 
cases give estimates for the activation barrier of 1.6 and 1.1 
eV, respectively. The experimentally measured barrier is 1.24 
± 0.08 eV.-^ 

on the connecting silicon atoms). 

By comparing the three panels of Fig. [S]we can now un- 
derstand the orbital origin of the ARPES bands. Band 
1' is the simplest and corresponds to the antibonding 
Da band. Band 2 is a superposition of bands S and Db, 
which are energetically for energies below the Fermi level. 
Band 1" has more complicated origin. It arises from two 
effects: strong rehybridization (related to the dimeriza- 
tion) around A2 of bands S and Db , and the energy shift 
discussed above from the adatom potential. 

We can also now better understand the two gaps cre- 
ated in Fig. [5] by adatom doping. The orbital origin of 
the those bands can now be assigned by noting, first, that 
the bands in Fig. OJa) are identical to those in Fig.[B]Jc), 
folded twice along the labeled vertical lines. It should 
be clear that the electronically induced gap illustrated in 
Fig. [SJb) arises from hybridization of bands S and Db, 
which cross near F when folded into the 5x4 zone. Like- 
wise, the adatom-induced gap illustrated in Fig. [5jc) is 
created entirely within the Da band at the second A4 
point, which folds into the A4 zone boundary of the 5x4 
zone. 



VII. DIFFUSION OF SILICON ADATOMS 

The silicon adatoms decorating Si(lll)-(5x2)-Au are 
not immobile.— Sequences of STM images of sur- 
faces prepared with adatom coverage near the equilib- 
rium value 1/8 show evidence for diffusion of adatoms 
along the Au-Si chain direction, by a defect-mediated 



mechanism of unknown origin 42. From measurements 
of the mean-square displacement at different temper- 
atures an effective activation barrier of 1.24±0.08 eV 
was extracted.*- In this section we examine whether the 
structural model proposed in Fig. [1] can account for this 
diffusion barrier. 

Although it is usually straightforward to calculate a 
barrier for the diffusion of an atom on a well-defined 
surface, Si(lll)-(5x2)-Au presents an interesting com- 
plication. Normally one studies the diffusion of a single 
atom in a supercell representing the clean surface, and 
relaxation of the surface around the diffusing atom is al- 
lowed. But for Si(lll)-(5x2)-Au the dimerization of the 
substrate is determined by the location (and coverage) 
of the adatoms whose very motion is under study. Thus, 
if the supercell is taken to be 5x4 and the substrate is 
allowed to relax, then the phase of the dimerization will 
track the position of the adatom. This "conforming sub- 
strate" is physically unrealistic and will underestimate 
the true diffusion barrier. 

A more plausible scenario, in which the phase of the 
dimerization remains fixed while the adatom diffuses, is 
difficult to implement without constraints or a larger su- 
percell. We take a simpler approach and adopt a per- 
fectly rigid 5x4 substrate, which will overestimate the 
true barrier. The true barrier must then be between those 
of the conforming and the rigid substrate. 

These two potential energy surfaces were calculated 
using a 5x4 unit cell in which the projected position of 
the adatom along the chain direction was constrained. 
For the conforming scenario, all other degrees of freedom 
were relaxed. For the rigid scenario, the two remaining 
adatom degrees of freedom were relaxed while all other 
atoms were were fixed at their original 5x4 positions. 

Figure [7] shows the two resulting potential-energy sur- 
faces. For the conforming substrate the potential-energy 
surface has 5x2 periodicity and an activation barrier of 
1.1 eV, while for the rigid substrate the periodicity is 5x4 
and the activation barrier is 1.6 eV. These two barriers 
nicely bracket the experimental barrier of 1.24 eV. More 
detailed studies using realistic boundary conditions will 
likely play an important role in unraveling the nature of 
adatom diffusion on this surface. 



VIII. OUTLOOK 

The structural model proposed here for Si(lll)-(5x2)- 
Au resolves one of the longest-standing unsolved recon- 
structions of the silicon surface. Predictions based on this 
model — for the dimerization of the underlying substrate, 
for the saturation coverage of silicon adatoms, for phase 
separation into adatom-covered and empty regions, for 
detailed STM imagery, for electronic band structure, and 
for the diffusion of adatoms — are in excellent agreement 
with experiments. Moreover, the physical mechanisms 
underlying many widely studied phenomena in this sys- 
tem have now been elucidated. 
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A number of other issues awaiting theoretical study can 
now be directly addressed. These include, for example, 
the suggestion that phase separation is accompanied by 
charge separation and the formation of a Schottky bar- 
rier at the interface between adatom-covered and undec- 
orated regions^ the structure and motion of the domain 
walls, within the Au-Si rows, that form when the spacing 
between two neighboring adatoms is not commensurate 
with the 5x2 substrate and the exploration of the fun- 
damental limits on using Si(lll)-(5x2)-Au to store and 
manipulate digital information at densities comparable 
to that of DNA4ii^ 

Finally, we anticipate a renewal of theoretical inter- 
est in Si(lll)-(5x2)-Au as a physical realization of a 
nearly one-dimensional metal. Most studies to date have 
used parametrized models because a definitive structural 



model has not been availablciSi With a model in hand, 
the door is now open to insights from new theoretical 
investigations as well. 
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